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Abstract 


The flowing-afterglow technique in conjunction with computer modeling of the flowing plasma has been used to 
determine accurate dissociative-recombination rate coefficients a for the ions O^, HCO + , CHJ, C H + , H 3 0 + CO + 
HC0 2 +, HN 2 0\ and N 2 0+ at 295 K. We find that the simple form of data analysis that was employed in earlier 
experiments was adequate and we largely confirm earlier results. In the case of HCO + ions, published coefficients range 
from 1.1 X 10 1 * * * * * 7 to 2.8 X 10 7 cm 3 /s, while our measurements give a value of 1.9 X 10“ 7 cm 3 / s - 


1. Introduction 

Rate coefficients for dissociative recombination 

(DR) of positive molecular ions with electrons are 

often needed to construct models of cold, weakly 

ionized plasmas, e.g. terrestrial and planetary iono- 

spheres [ 1 ] and interstellar gas clouds [2,3]. Many 

ions have been studied subsequent to the work of 
Biondi at al. who used the stationary afterglow tech- 
nique [4-10]. More recently other techniques, espe- 
cially flowing afterglows [11-18], but also merged 
beams [19,20], and ion storage rings [ 21 ] have be- 
come important tools for such work. The chemical 
versatility of flowing afterglows has made this tech- 
nique particularly productive, but a rigorous analysis 
of the data obtained by this method is a rather 
difficult task. In this article we present measurements 
of the recombination coefficients for several poly- 
atomic ions obtained using a flowing afterglow/ 
Langmuir probe (FALP) apparatus. In part, the pre- 
sent study was motivated by the recent finding [15,22] 
that the de-ionization of H 3 and D 3 + afterglow 


plasmas deviates significantly from the simple be- 
havior that is expected for binary recombination of a 
plasma containing a single ion species. The com- 
monly used approximate methods of deducing 
recombination rates from graphs of the reciprocal 
electron density versus time or position tend to mask 
deviations from a simple recombination decay, and 
one can never be sure that the inferred recombination 
coefficients are truly constant during the afterglow. 
We have sought to improve this situation by devel- 
oping better approximations for the data analysis 
which would reveal subtle deviations from a simple 
recombination decay. 


2. Experimental 

Our flowing afterglow technique [22] and the use 
of Langmuir probes for electron density measure- 
ments [23] have been described in detail elsewhere 
and will only be summarized here. A microwave 
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discharge in helium is used to generate helium 
metastable atoms. He \ which further downstream 
ionize argon by Penning ionization. In this work, 
several of the desired ion species were formed by 
charge transfer from Ar“\ 

Ar + + XY — > XY + + Ar , (1) 

where XY is 0 2 . C0 2 , or N 2 0. To make C 2 H7 
ions, CH 4 was added to the afterglow leading to the 
reactions 

Ar 4 + CH 4 — >CH.| + H + Ar (2) 

CHt + CH 4 -► C 2 Hj + H 2 . (3) 

Protonated molecules were produced by first making 
H ; and subsequently transferring the proton to the 
molecule under study, i.e., 

Ar + + H 2 -» ArH + + H (4) 

ArH ^ + H 2 — ► H t + Ar (5) 

H; + XY -> HXY 4 T H 2 (6) 

where XY is CH 4 , C0 2 , N 2 0, CO, and H 2 0. The 
quadrupole mass spectrometer at the downstream 
end of the tube was used to monitor the ion composi- 
tion of the plasma and to ascertain that no significant 
impurity ions were present. 

The movable Langmuir probe measured the elec- 
tron density as a function of distance ~. Under 
typical experimental conditions the electron densities 
in the upstream region of the flow tube were of the 
order of 5 X 10 10 cm" 3 , the flow tube pressure was 
0.9 Torr. and the central velocity was ~ 5500 cm/s. 

2.1 . Data analysis 

The usual practice of deducing recombination co- 
efficients from the slopes of the inverse electron 
density versus distance z plots is justified only if 
ambipolar diffusion in radial and axial directions can 
be safely ignored. The on-axis electron density then 
follows approximately the equation 

— — = ! + «-, (7) 

*c(0 n«U = °) v 

where v is an ‘effective flow velocity’ whose precise 
relation to the central How velocity does not seem to 
be well defined. While it is true that the initial slope 
of a graph of Eq. (7) using experimental values for 


n c often yields good values for the recombination 
coefficient, the problem is that further downstream 
the transition from a recombination-controlled to a 
diffusion-controlled afterglow leads to a gradual in- 
crease in the slope. In the absence of a better method 
of analysis, it is difficult to distinguish this effect 
from an actual change of the apparent recombination 
coefficient which may be present. Such changes may 
result from changes in the ion composition, quench- 
ing of vibrationally excited ions, or contributions of 
electron-stabilized recombination. 

We have developed an accurate numerical model 
of the afterglow processes [22] that includes the ion 
chemistry and recombination, mixing of reagent 
gases, diffusion and transport of all species in both 
axial and radial directions in a flow with a parabolic 
velocity profile. In practice, however, the calculation 
time was found to be extremely long for routine 
deduction of recombination coefficients, and the 
model was used mainly to validate a simpler, approx- 
imate method of analysis. 

In the approximate solution of the electron conti- 
nuity equation the gas flow velocity is taken to be 
radially independent and the radial variation of the 
electron density is approximated by the fundamental 
diffusion mode in a cylinder. With these simplifica- 
tions the continuity equation for the on axis electron 
density is of the form 


d« e ( r = 0) 


1 


a 2 n e 

n \ 


-an] T D a | 


L 

"o 

a; 2 

*J] 


(«) 


where A = /?/2.405 is the fundamental diffusion 
length for a cylinder of diameter /?, D a is the 
ambipolar diffusion coefficient and u 0 is the central 
flow velocity. A fourth-order polynomial least- 
squares fit to the measured c-dependence of 1 / /? e is 
inserted into Eq. (8), the needed derivatives are 
obtained analytically, and the equation is then solved 
for a(c). If the resulting a(z) is independent of z, 
then it can be concluded that a single, constant 
recombination coefficient prevails during the after- 
glow. This approximate method of data analysis was 
tested by creating computer-simulated data using the 
accurate numerical model and then analyzing them 
using the approximate method. The results confirmed 



7. (iougousi et at. / Chemical Physics Letters 265 1/997) 399-405 


401 


that the approximate method is entirely satisfactory 
for analyzing experimental data. 

Ambipolar diffusion coefficients for the different 
species are needed as input parameters for Eq. (8). 
For most ions they were calculated from the known 
mobilities [24-26] in argon and helium using the 
Einstein formula relating the diffusion coefficients to 
the mobilities. In cases where the mobilities in argon 
were not known the Langevin formula [27] was used 
to obtain an estimate. The central velocity v 0 was 
calculated from the known gas flow rate and the total 
pressure. 


3. Results and discussion 

The recombination coefficients were measured for 
the ion species listed in Table 1. Before data were 
recorded, the mass spectrum was taken to verify that 
the desired ion species indeed dominated the plasma. 
As a test case for the technique and the data analysis 
the recombination coefficient of Ot ions was mea- 
sured. A value of ( 1.9 ± 0.2) X 1()“ 7 cm 3 /s was 
obtained, in good agreement with the generally ac- 
cepted value of (2.0 ± 0.2) X 10 7 cm 3 /s. Fig. 1 
shows an example of the electron-density decay curve 


Table I 

Dissociative recombination coefficients tor polyatomic ions measured using a flowing afterglow/Langmuir probe apparatus. Coefficients 
measured in earlier experiments are also given 


Species 

Present work a (10 7 cm 3 /s) 

Previously reported values 
or (10 7 cm 3 /s) 

Reference 

Method a 

o; 

1 .9 ± 0.2 

1 .95 

Mehr [4] 

SA 



2.20 

Kasner [5] 

SA 



2.30 

Mul [19] 

MB 



1.95 

Alge [11] 

FALP 



1.90 

Walls [28] 

IT 



2.00 

Spanel [ 1 3] 

FALP 

MCO 

1 .9 ± 0.3 

2.0 

Leu [8] 

SA 



1.1 

Adams [17] 

FALP 



2.4 

Ganguli [10] 

SA 



2.8 

Amano [29] 

IR 



2.2 

Rowe [ 1 4] 

FALP 



1.5 

Smith [ 15] 

FALP 

co; 

3.5 + 0.5 

3.8 

Weller [6] 

SA 



4.0 

Gutcheck [30] 

SA 



3.1 

Geoghegan [ 1 6] 

FALP 

n.o 1 

3.6 ± 0.6 




hco; 

4.6 ± 0.7 

3.4 

Adams [12] 

FALP 

hn 2 ct 

4.5 ± 0.7 

4.2 

Adams [12] 

FALP 

H ,O k 

8.0± 1.5 

10.0 

Herd [18] 

FALP 



3.4 

Mul [20] 

MB 



13.0 

Leu [9] 

SA 



20.0 

Heppner [3 1 ] 

IT 

ch; 

9.0 ± 2.0 

11.0 

Adams [ 1 7] 

FALP 



14.0 

Smith [ 1 5] 

FALP 

c 2 h; 

9.0 ± 2.0 

7.4 

Adams [ 1 2] 

FALP 


J SA: Stationary afterglow. FALP: Flowing afterglow/Langmuir probe. 1R Infrared absorption. IT: Ion trap. MB: Merged beam. 



402 


T. Gougousi et al. / Chemical Physics Letters 265 (1997 ) 399-403 



0 2 4 6 8 10 12 14 

distance (cm) 


Fig. 1. Values of a(z) (Filled squares, right-hand ordinate) from 
measured n c ( z) data (open circles, left-hand ordinate) for an 0 2 + 
plasma. The solid line through the n e (z) data represents a fourth- 
order polynomial Fit. The solid line through the cr(^) values 
represents the average, while the dashed lines indicate the one 
standard deviation interval. 

and the deduced a(z) for ions. As was ex- 
pected, a(z) remains constant over a factor of ten in 
electron density, supporting the conclusion that 
deionization occurs by binary recombination of a 
single ion species. Data of equal quality (see exam- 
ples in Fig. 2 and Fig. 3) were obtained for most of 
the other ions in Table 1. As a rule, the a(z ) values 
for all ions were constant within 10 to 15% from the 
mean. None of the ions exhibited the strongly 
anomalous decay that has been observed for H 3 ions 
[15,22]. 

For most ion species our recombination coeffi- 
cients agree weTT with earlier data where such are 
available. The agreement among various determina- 
tions of a(HCO + ) is poorer than for other ions. The 
lowest value was reported by Adams et al. [17] at 
1.1 X 10 “ 7 cm 3 /s and was revised later by Smith et 
al. [15] to 1.5 X 10 “ 7 cm 3 /s, while the highest value 
(2.8 X 10 “ 7 cm 3 /s) was reported by Amano [29]. 
With the exception of these two cases the rest of the 
experiments yielded values near 2.2 X 10 “ 7 cm 3 /s 
and our result 1.9 X 10 “ 7 cm 3 /s is in good agree- 
ment with them. 

The state of internal excitation of the various ions 
is generally not well known. Our data contain no hint 
that the recombining ion species consisted of mix- 



0 2 4 6 8 10 


distance (cm) 

Fig. 2. Values of a(z) (Filled squares, right-hand ordinate) from 
measured n c O) data (open circles, left-hand ordinate) for an COT 
plasma. The solid line through the n t {z) data represents a fourth- 
order polynomial Fit. The solid line through the <*(?) values 
represents the average, while the dashed lines indicate the one 
standard deviation interval. 

tures of different vibrational states and that these 
recombined with significantly different recombina- 
tion coefficients. If this were the case, the finding of 
a constant recombination coefficient would be rather 



0 2 4 6 8 10 12 14 16 

distonce (errs) 

Fig. 3. Values of ar(z) (Filled squares, right hand side axis) from 
measured n e (z) data (open circles, left-hand ordinate) for an 
HCO + plasma. The solid line through the n c (z) data represents a 
fourth-order polynomial Fit. The solid line through the cr( z) 
values represents the average, while the dashed lines indicate the 
one standard deviation interval. 
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surprising. Since vibrationally excited ions are often 
quenched to lower states in collisions with molecular 
gases, the presence of excited ions might also lead to 
an apparent dependence of the recombination coeffi- 
cient on reagent pressure. However, measurements at 
different reagent concentrations resulted in no more 
than 10 to 15% variation of the recombination coef- 
ficient indicating that either the ions were in their 
ground state or that there is little dependence of a 
on vibrational state. 

While our data agree well with many other deter- 
minations, it should be noted that most other deter- 
minations were made in plasmas of similar pressures 
and electron density. There is still lack of accurate 
absolute recombination coefficients (or cross sec- 
tions) measured under essentially single-collision 
conditions, such as in merged-beam or ion storage 
rings. 
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